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Abstract 

The rate for the final state radiation of hadrons in leptonic Z decays is evaluated, using as 
input experimental data for a{€^e~ — > hadrons) in the low energy region. Configurations 
0> ' with a lepton pair of large and a hadronic state of low invariant mass are dominant. A 

relative rate ^ilhad/^ii — 6-3 x 10~^ is calculated. This result is about twice the prediction 
Qh' based on a parton model calculation with a quark mass of 300 MeV. The rate for secondary 

production of heavy quarks is calculated in the same formalism. 
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00 . Introduction. Rare decay modes of the Z boson into final states with two fermion-antifermion 

^ pairs could provide the first signal for new physics, in particular for the production of a new 

O . scalar boson. They are, however, also predicted in the context of the Standard Model. A 

Q^ ! thorough understanding of production rates and distributions is, therefore, mandatory. The 

"^ ! final states can be sorted into three classes: Into a class (i) with quarks only, another class (ii) 

Ph! which involves a quark plus a lepton pair and finally, a class (iii) with leptons only. Class (i) 

P-i contributes in order a^ to the purely hadronic final state and is part of the hadronic decay rate 

r-| ; calculated in |I], g, ^. Class (iii) can be treated by standard methods of perturbation theory. 

>■ ' These formulae have also been applied to P class (ii) final states, assigning a mass of 300 MeV 

KA to the light quarks. However, in practice the dominant contributions consist of a few pions with 

low invariant masses. A treatment, based on the experimental knowledge of Rhad-, the rate for 

hadron production through the virtual photon, should complement and in fact superseed the 

purely perturbative treatment. 

The techniques employed in this paper have been developed |Q in the context of initial state 

radiation. They are based on the observation that the information relevant for the production 

of low mass hadronic states through the virtual photon can be encoded in a few moments 

of Rhad- For hadrons these have to be calculated numerically. For leptons, however, they 

can be calculated analytically, providing a convenient test of results that can be found in the 

literature. 

The subsequent treatment will be limited to contributions induced by virtual photons only. 

Contributions to the rate where the secondary fermions are produced through neutral current 

interactions can safely be neglected. 
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Figure 1: Typical diagram, describing the radiation of hadronic states from the II final state. 



Hadronic final state radiation. Neglecting the mass of leptons, the rate for Z — > // + hadrons 
which is induced through the amphtude depicted in Fig.l can be cast into the following form: 
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where Li2 denotes the dilogarithm Li2(y) = — Jq '^ ,~^M x. Neglecting trivial factors, F(s'/s) 
gives the rate for the decay of a vector boson of mass s into a vector boson of mass s' plus 
a pair of massless fermions (// ) . The formula is valid, as long as the fermion mass mi is far 
smaller than vs'- Eq.|l] can be viewed as a representation of the rate for II + hadrons through 
a superposition of a series of resonances with masses s', distributed with weight R{s'). 
Since the threshold for pion production is located at 2mT^ and the minimal s' is far larger than 
Amf this approach is well justified if / is an electron. The approximation can also be applied for 

which is significantly larger 
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muons, since R{s') starts to contribute significantly for ys 
than the muon mass. However, it is not applicable for t^t^ + hadrons. It is valid for the 
contribution to the rate for I1I1I2I2 where the heavier lepton is radiated off the light leptons, 
and more generally, if the invariant mass of the " secondary" lepton pair is constrained through 
appropriate cuts. 

The final state under consideration consists typically of a lepton pair with large and a hadronic 
state with low invariant mass which is emitted collinear with the lepton momentum. [The 
converse possibility, where the virtual photon is radiated off the quark can be formally treated 
with a similar technique. However, the final state consists mostly of relatively soft and collinear 
photons of low invariant mass and the corresponding hadronic amplitudes are difficult to tackle. 
Since this configuration leads, typically, to a lepton pair of low invariant mass, it can be easily 
distinguished from the former configuration.] 



The evaluation of (|I]) proceeds as follows. The integral can be split into a part which is governed 
by the large s' behaviour of R and a remainder which is sensitive to the details of the threshold 
behaviour: 

r ^R{s')F (s'/s) = R{oo) t ^ , ^F{X) + T^^ [R{s') - i?(oo)] F (s'/s) , (3) 
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with 



R{s) '-=^ R{oo) + C(ln-^ s) . (4) 

At this point it is useful to observe that for /i one may take any positive value, as long as the 
function R{s) vanishes by definition below 2m,r • However, for the following discussion /z will 
be simply identified with m^ . The first integral can be evaluated in a straightforward way. 
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where the trilogarithm Lis is defined by Liad/) = J^ ^^^dx, (2 = Li2(l) = ^, Cs = Li3(l) 

1.2020569... . 

For large y = s/(2yu)^ one obtains 
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The second integral is dominated by low s'. The function R{s') approaches rapidly the asymp- 
totic value R{oo) and the remaining integral can be extended to infinity. Only small values of 
s'/s contribute to the integral, and F{X) can be approximated by its behaviour for small A 
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Combining the two contributions one arrives at 
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where the moments Rn are defined through 
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This result can be used to predict the inclusive hadronic rate, if we employ R{s) as provided from 
experiment p. A convenient approach adopted to the experimental procedures is to decompose 
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Figure 2: Parametrisation of R{s) in the threshold region, without resonances, as provided by 

/!/■ 



Table 1: Input parameters used for the narrow resonances and relative rates of hadronic radia- 
tion. 



Channel 


Fe+e- [keV] 


M [MeV] 


Relative Rate 


2n 






2.55 X 10"^ 


Continuum 






2.86 X 10"^ 


UJ 


0.60 


781.95 


2.34 X 10-5 





1.37 


1019.41 


3.56 X 10-5 


^ 


5.36 


3096.93 


2.25 X 10-5 


f 


2.14 


3686.00 


6.55 X 10-6 


T 


1.34 


9460.32 


5.85 X 10-^ 



the final states into continuum (Fig. 2) and resonance contributions. For the continuum one 
finds 

R{oo) = 11/3 , Ro = -13.240 , Ri = 20.875 , R2 = -3.587 . (11) 

For the two pion (=p) contribution the parametrisation derived in P| (set 1) is employed which 
leads to the following moments: 



R{oo) = , Ro = 4.319 , i?i = -8.422 , R2 = 8.547 . 
For the remaining resonances the narrow width approximation can be applied 
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The input parameters are listed in table |1|, together with the predicted relative rate. In total 



Table 2: Contributions of different energy ranges to the relative rate ^uhad/^ii- 



Energy range 


Contribution to r^^-^^^^/F^^ 


< ruhad < 2 GeV 


3.88 X 10-^ 


2 GeV < nihad < 3 GeV 


5.87 X 10-5 


3 GeV < TTihad < 10 GeV 


1.52 X 10"^ 


10 GeV < TTlhad < Mz 


3.16 X 10-5 



we find the following result 
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It should be reiterated that this contribution corresponds to configurations where the invariant 
mass of the lepton pair is typically large and the mass of the hadronic system small. This is 
evident from table |^ where the contribution for nihad > 10 GeV is also shown separately. 
The numbers listed in tables |I| and are difficult to compare with those from 0] . There the 
calculation was performed numerically for a qqll final state with a fixed quark mass of 300 MeV. 
For / = e and to lesser extend for / = /x their result is necessarily dominated by final states 
with large hadronic mass, that is a multi-hadron final state plus a lepton pair of low invariant 
mass. It seems experimentally difficult to isolate these events among the multi-hadron states, 
which occur with a relative rate below 10"'^. It is however an interesting exercise to apply (^ 
and the formulas discussed below to llqq production where again the hadrons are radiated off 
the lepton. These can experimentally isolated by considering events with large invariant mass 
of the lepton pair. Adopting the quark mass values of |@, viz. m„ = rud = rus = 300 MeV, 
rric = 1.5 GeV and m^ = 4.5 GeV one predicts 
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5.259 X 10" 
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2.858 X 10 
3.938 X 10"^ 



2.256 X 10" 



(15) 



We find that this model which has been used in 0, underestimates the rate for radiation 
by a factor 0.63. An experimental analysis based on this model could, therefore, lead to 
drastically wrong conclusions. The reason for the discrepancy is fairly evident: The events are 
dominated by hadronic states with low invariant mass, where the parton model provides a poor 
approximation. This is reflected in a strong sensitivity of the results towards the choice for m^. 



Four lepton final states. The same formalism can also be applied to purely leptonic final states, 
using the analytical results for the moments^ H] 
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The numerical results for T^+^-ii/Te+e- are 3.319 x 10 ^ for / = /i and 3.214 x 10 ^ for / 



T. 



It should be reiterated that this result refers to final states where the // is radiated off the 
e"'"e~ pair and hence is dominated by // with low invariant mass. Eq.^ can be also used to 
predict the rate for lepton pair radiation if a cut Si+i- > sq ^ Amf is employed. In this 
case one may simply put Rq = Ri = R2 = and y = s/sq . Characteristical examples are 
T^+^-i+i-/Te+e- = 5.321 X 10-71.268 X 10-78. 716 x 10"^ for ^ = 0.1/1/10 GeV respec- 
tively. 



Secondary radiation of heavy quarks. Replacing the factor (a/vr)^ by (as/7r)^2/3, the pertur- 
bative result of eq.^ combined with eq.|l^ can be applied to predict the rate for secondary pro- 
duction of heavy quarks in Z — > gg events, where q denotes a light quark. For mi, = 4.5/4.8/5.0 
GeV one predicts a relative rate of 1.545 x 10-^/1.391 x 10-'^/1.301 x 10-'^. The corresponding 
numbers for secondary charm production [nic = 1.3/1.5/1.8 GeV) are significantly larger and 
amount to 8.339 x 10-77-107 x 10-75.729 x 10-^ respectively. These rates should be well 
observable with present day statistics. 



Combination with virtual corrections. The third and second power of the logarithm of y cancel 
after combining real radiation of hadrons (or leptons) with the corresponding virtual corrections 
to the form factor H 
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and one obtains for the correction to the total rate 
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The coefficient of the remaining logarithm can be easily understood: The lowest order correction 
(from real and virtual photon radiation) is given by a factor (1+|^). The fine structure constant 
a can be replaced by the running a(s) = a/(l — n(s)). The vacuum polarisation from hadronic 
(or leptonic) intermediate states can be expressed in terms of R through 
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as f°° ds' R{s') 



IT 3 J4/i2 s' S — S' 



t The result for the total cross section in leading logarithmic approximation can be found in 



and one obtains for the complete correction (without 0{a'^) photonic terms) 

3a(s) 1 /a\2 / 23\ 

i + i^nb) r'-TJ- P") 

An amusing test of ([T8|), combined with the moments from (p!6D, is provided by the results of 
|I|]. If one replaces the SU{3) factors in this paper by the corresponding abelian coefficients, 
and, furthermore, relates the MS coupling aj^ to the conventionally defined fine structure 
constant 

a = a{l + al\n^] (21) 

\^ 6 m'^ ) 

one obtains agreement between the two results.| 
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